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SUMMARY 


The heat flow measurement system obtains a temperature profile 
trom the surface of the test structure by means of thermocouples. The 
thermocouple signals are used as inputs to an analogue computer which 
solves the appropriate heat flow equation. This computer is time shared 
between the 4 channels of heat control by means of a high speed stepping 
Switch. The system is very versatile since it is independent of the 
source of heat. Thus it is inherently more exact than a radiation trans- 
ducer system because it accounts for convective and other losses which 
occur in the simulation process. It has the added advantage that it can 
be used to measure cooling rates when the cooling results from con- 


vection as well as from radiation. 
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I. INTRODUCTION 


The structures laboratory of the Department of Aeronautics and 
Astronautics at Stanford University concentrates its research in the field 
of thermal effects. Thus a Structural test facility which includes the 
ability to Simulate aerodynamic heating is of prime importance. The 
basic deterrent to a university owning and operating such a facility is 
that, if the facility is of reasonable size, the cost is normally prohibi- 
tive. The problem wes therefore reviewed with the intention to design 
and construct a facility that would be economical enough for a university 
to own and reliable enough to get good research results. It was desired 
that the facility should have a capacity of a megawatt or more, and that 
it should be automatic in operation. It was considered essential that it 
should be computer controlled, and that the loop closing should be based 
on a measure of heat flow and not on temperature or some other property. 

For a control system of this kind, it has been said that, "....(the) 
weak link of all methods at the present time is the assessment of the heat 


input. The fact that cifferent methods are adopted for this purpose when , 

different forms of heating are used complicates the matter still further. nt 
How this weak link has been strengthened in the construction of a 

versatile heat flow measuring system compatible with most forms of 


heating is the subject of this thesis. 


*Superscripts refer to the reference list appended to this report. 


z 





Pee ene 2 Ol SY stig 


The complete control system is functionally diagrammed in 
Figure 1. Simulation of the convective heat transfer from the air to 
the unit under test is achieved by a solution of the Newtonian equation 


ior cooling. : In symbols this equation is 
a(t te 
Ga My  tgut 4 ote 
where 


q = heat flux/unit area (BTU/ft7br) 
h = the heat transter coefiicient (BTU/ft"br°F) 


ae ne adiabatic wall temperature (o>) 


O 
te = Surface temperature (F) 


- total radiation pickup (BTU/ft“br) 


q.. 


4. 
ot. - Stefan-Boltzmann loss (BTU/tt“hr) 


A digital program ofh, t__., and q_ (external radiation from 


WwW 
the sun and nuclear explosion) is derived from a proposed flight profile. 
This program is converted to voltage signals by the digital to analogue 
MCONVYCItCL meine VOlLage Ley is combined with a voltage analogous to 
surface temperature in»a summer and the result multiplied by h ES 

This product is added to q.. to give the desired heat flow signal, q qd: 
From the feedback loop the measured heat flow to the test specimen, 

da? is contrasted to gq 4 in acomparator. The output of the comparator, 
c , causes the saturable core reactors to vary the power supplied to 

the radiant heaters. The control loop is closed by the heat flow measure-~ 


oe : cine 
ment system” and the correct heat flow is maintained. 








Automatic operation of the facility is possible trom the point 
at which a flight profile, procrammed into digital form, is fed into 
a digital to analogue converter. An analogue computer receives 
the analogue signal from the converter and computes a command 
signal to adjust the power supplied to the heaters. The feedback 
systems then produces a signal proportional to the radiated heat 
absorbed by the test specimen. This signal is used to close the 


control loop and render the system automatic. 








1 Teo ee aL ING TEE SYSTEM OPERATION 


The Fourier equation for conductive heat flow within a body, where no 
convective heat flow is possible, can be combined with the principle of 
; 5 ; ee? . 
conservation of energy. This combination for a thin block at a surface 


(Figure 2) subjected to a heat flux, q, is outlined below: 


in the x direction 


heat flow in = q(dy dz)dt (1) 

heat flow out = -K(dy dz)(6/6 x)(T + dT/6x~- dx)dt (2) 
where 

T = temperature 

t =time 

K = material conductivity 


in the v cirection 
heat flow in =-K dz dz(dT/dy~- di) (3) 


heat flow out = -K dz dx(d/dy)(T + dT/dy - dy)dt (4) 


in the z direction 


heat flow in = -K dx dy(8T/dx - dt) (5) 

heat flow out = -K dx dy(d/0z)(T + OT/3z - dz)dt (6) 

the heat stored = cp dx dy dz(dT/odt - dt) (7) 
where 


po = material density 


c = specific capacity of the material 


Combining equations for tne heat in = the heat out + the heat stored 


and then cancelling like terms results in a simpler equation, namely 





x 2 


» a ») 
Na 
pe ST dx=-K oT -K aT at P+ So | dx +q (8) 
Ox a: OZ 


For a like block below the surface the equation is 


. 2 2 2 
po ST ax = K(9T 4 ST ST) ay (9) 
ails ox Sy 32 


Summing the set of equations results in 


nal 
pec { 3 2 dx = G «(2 Ag. dx+ K ST +Q 


O O Ss oy Aig ~ x=0 


(10) 
wnere 


n = numoer of divisions 


I 


spacing 


This equation holds for a slab of material whicn is bounded by the upper 
surface of the structure and a plane parallel to that surface so long as 
this bound is not the lower surface of the skin. This restriction is ne- 
cessary to avoid the undefinable loss from the inner surface due to radiation 


and convection. 


Equation (10) can be re-arranged for an expression for heat flux. 





na na 9 9 A 
a= pe | ar _ { «(2S 2G +See 
O (© O Ox ay OZ ox =e. 
ha na 9 
= pe | aT ax + | (2. ee ie Ok (11) 
CO enon O s 2 OX x=na 
i 3z 


The spatial derivatives (3°T/3y" oe ie 20/32") at any level in the 


2 2 e } ° 
direction of the inward normal can be expressed in the form 








Q/ 


mn I 
a 1 =¢ 
2 C (aT, pT, Hee 77] (12) 
where C,@, 6, y .2ssare cocificients as reproduced in Table l 


from Reference 6, 


The integrations to be made can be done by an approximation. By 


application of Simpson's rule for n= 2 equation (12) becomes 


2 


Z 
y 


G7 
a 


1 
Se [Eo ene ry ae (13) 


Q/ 


where the equally spaced temperatures are measured in the y-z plane 


as depicted in Figure 2. 


Equation (11) then can be written 





f Gy em ee Gs ae ee 
q ii o ena St Fa 30 ( on a 7 9) or 39 ( O i a ats 9) 
oT; | 
= C vas Colo +) Eo ee i to) (14) 


All the terms in equation (14) are directly determinable from a spatial 
temperature survey. The differentiation can be performed in an analogue 


computer as can the summations and multiplications. 


The preceeding analysis outlines the method of the heat flow solution 


utilized in the feedback computer system for all 4 channels of control. 





Ve OE eo ON eect GLa NN CL 


Eighteen thermocouples are capacitance-welded to the surface 
of the test structure. The temperature survey used in the first channel 
of control appears in Figure 3. The voltages generatec by the thermo- 
couples are connected to a hieh-speed stepping switch which sequentially 


connects the voltages to the appropriate computer elements. 


The analogoue computer consists of a d-c power supply, 4 chopper- 
stabilized amplifiers, 1 differentiator, 1 summer, and 5 peak follower 


amplifiers. Details of these components are discussed in a later section. 


To clarify the different functions of the computer, a step-by-step 
Solution of the meating equation for channel 1 is pres@nted (s@e Ficure 4). 
T, , aS a voltage, is received on the first step by bank number 1. From 
bank 1, T, is amplified and then age a to bank 2. From here T, 
is again amplified and then dif fereniated. TRE Ourour ot Lhe Pe ae 
is temporarily stored in capacitor, C. . This charge is drawn from Cy 
on the second step of bank 3 by peak follower 1 and then connected to 

a eo : : 


=e -T, are 
joined to banks 4, 5, 6, 7 and 8 resnectively. These signals are 


summer S. Also onthe second sten T a% } 
summed and amplified in chopper amplifier 1. Through bank 9, the sum 
goes to chopner-amplifier 2, ater amplification this voltage is diverted 
through the second step of switch 2 tosummer §. This summer sums 
the voltages as required by equation (14) . This solution, q, is 
switched through banks 10 and 11 to storage in capacitors Co and Cy. 
Two banks are necessary here so that the total q may be decreased as 
well as increased. This is done by grounding Cy through bank 11 on the 
third step. Bank 12, still on the second step, diverts q through peak 


follower 2 to the comparator thus completing the control loop. 


It should be noted that the radiation losses from the test specimen 
have yet to be taken into account. To compensate for this effect, it was 


4A 
determined that a oT) term, where o is the,Stephan-Boltzmann 





constant, could be apnroximated by $ straight line slopes. (For this 
approximation see Figure 5.) Accordingly, a biased diode circuit 
was devised to follow the approximation. (See Figure 6) T; is first 
amplified to a workable voltage and is then summed with a ee... 
voltage to make it propvortional to absolute degrees Fahrenheit since o 
involves absolute temperatures. The output from the summer is then 
acrnitted at the proper level by the biased diode circuit to achieve the 
7 Changes. slope necessary to approximate the fourth power function. 
It is to be noted that at smaller voltages, i.e. lower temperatures, the 
radiation compensating circuit is opened. At low temperatures, bow= 
ever; the radiation loss is negligible (Figure 5) and has virtually no 
effect upon the total equation. The amplified voltage output from the 
radiation compensation circuit is attenuated at summer S (Figure 4) 


to the proper level to correct q. 








Vow colours o DAITS 


The thermocouples used were made of 28 gauge chromel-alumel 
wire junctioned with a bead weld? Standard chromel-alumel thermo- 
couples were found to give good response up to temperatures of 2450° 
Fahrenheit. Although platinum-rhodium thermocouples extended the 
upper temperature limit to 3100°F , it was felt that the additional 
expense involved was not worth the increase in temperature range. 2 
The thermocouples were capacitance welded to the test surface. The 
entire test area, including the thermocouples, was painted black in order 
to minimize the amount of radiant power required. Several commercial 
heat resistant paints were testec and the best found to be common Stove 
blacking. This paint was extremely resistant to heat. It was easily 
applied by brush and produced an even coat which unfortunately scratch- 


ed readily. 


A Philbrick compound regulated dual power supply was used to 
drive the computer amplifiers. This power supply has good stability 
and reliability. Its output of * 300 volts d-c is stable to * 100 
millivolts per day. ‘ 

The computer amplifiers used were either Philbrick K2X octal 
based plug-in amplifiers or a combination of the K2X andthe K2P. 

he K2X's response time is 1 microsecond and has a rated adrift of 
+ § millivolts per cay. Tue input impedance of each amplifier is 100 
megonms and the nominal output current is of the order of 4 milliamps. 
The amplification of the K2X is limited to 10 by stability criteria; 
consequently, K2X's were used aS summers and as peak followers. 
A peak follower is essentially an amplifier with unit feedback. It is 
used in conjunction with a temporary-storage capacitor and will provide 
an output voltage which is the peak value of the past history of the ca- 
pacitor charge. To attain higner stability when amplifying by 1000, a 
K2X amplifier is stabilized by a K2P chopper. rift rates for this 
setup are of the submillivolt level and accuracies in the order of 0.1% 


are attained. 





Toe canacitors used aS storace banks are all 0.25 microfarads 
With a2 time constant of 10 microseconds. This means that the capacitor 


is charged to 99% of its final value in that time. 


The device which enables the computer to time share its amplifiers 
is atyne 45 rotary switch manufactured by Automatic Electric. It is a 
heavy duty, raaenet-driven stepping switch with positive non-overthrow 
steoping. It can be operated by remote control pulses or by self-inter- 
ruption as is done in this application. The wipers, with wiping tips at 
both ends (180° Dpart),sonemmocd a im one Cirection over a semicircular 
arrangement or 12 banks of 25 contacts each. The first 3 banks are 
gold nlated ana are used in the differentiator circuit where the decrease 
in noise level is of prime importance. The switch is driven bya 24 
volt d-c power sundly to step at 56 times per second. Since the wiper 
is double ended, each contact is scanned at a rate of 112 times a second. 
Dwell tirce on a contact is 6.4 milliseconds or 80% of the time for each 
step. The switch is so wired as to complete the computations for one 
channel in 3 steps. The fourth step begins channel number 2 compu- 
tations and so on in sequence. Since the switch contains 25 steps, the 


channels are scanned 4 times for each revolution of the switch wipers. 


=e 





Vi. CONCLUDING REMARK 


The heat flow measuring system incorporated into the radiant 
heating structures testing facility has the primary advantage of being 
independent of the source of heat. It can be used to measure heat flux 
in systems utilizing convection, radio frequency induction, or radiant 
heating methods. It also can bo used to measuro heat flow in the cooling 
phase of a flight profile which could produce serious thermal transients 


comparable to those in the acceleration phase. 


Experiments with the heat flow measuring system described in this 
thesis were So successful that it was decided to extend the system of 
control to include the deceleration as weli as the acceleration phases. 

For this purpose, it is essential to introduce some form of cooling in 
addition to radiant cooling. After due deliberation, it was considered that 
the most successtul means of additional cooling would be to blow cold gas 
onto the surface of the specimen. I was decided that liquid nitrogen 
would provide the best heat sink; and as a conseauence, the heater arrays 
wers adapted to include cooler arrays (see Figure 7). These cooler arrays 
consist of a series of quartz tubes identical in character to the envelopes 
of the infra-red lamps. They were placed in a plane parallel to that of 
the neater laraps and had slots and holes along their lower surfaces. 
Toroughn these tubes was blown cold nitrogen from a gas generator inte- 
eral with the reflector. This arrangement did not interfere in any way 
with the passage of the radiant heat to the svecimen Since the quartz 
tuses cof the cooler array transmitted the radiation with no loss of ef- 
ficiency. The presence of the heat excnanger in association with the 
reflector surface gives two distinct advantages: first, the cold fluid is 
conveyed as near as possibie to the test specimen in the form of a liquid 
and second, the cooling of the reflector itself makes it possible to use 
intensities of infra-red not practical without this feature. In order that 
the systeza of regulation should not require additional servos for the 


I 
cooling case, overcooling was used and the balance of heat adjusted by 


ge 


means of tne intira-red system. fureacr cCevelopmental work alone these 
lines is clearly desirable, but it seems quite certain that the feasibility 
of such a method has been established. 

A second advantage, often oz the greatest importance, is that the 


system can be economically constructed from standard components. 


cial transducers are needed to measure the heat flow. 
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THE CALORIMETER AND TEST DATA 
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THE CALORIMETER AND TEST DATA 


For calibration of the transducers, a calorimeter was made from 
stainless steel. It consisted of a 1/32 inch thick disc, 3/4 inch in 
diameter. A diagram of the calibration equipment is given in Figure Al. 
The calorimeter was guarded along its edge by a ring of the same material 
and on its rear face by another disc. In this way the heat losses from the 
calorimeter were made negligibly small. This was verified by recording 
thermocouples attached to various points on the rear face of the calori- 
meter. Tne temperature of the disc was found to rise uniformly when 
subjected to a radiant field. A series of temperature versus time records 
was made for the calorimeter in the uniform field. (See Figure A2 for 
a picture of the test arrangement.) The tests were performed with the 
calorimeter at the same distance from the source that a test structure 
would normally be located. A radiation transducer was mounted in its 
normal position on the heater array used for the test. The output of the 
transducer was recorded simultaneously with the calorimeter temperature. 
The result of each test was thus a transducer output which corresponded 
to the calculated heat absorption rate. During actual structural test, the 
Same transducer signal will represent the same heat absorption rate by 
the test structure. One transducer was fully calibrated in this manner 


and the others calibrated against this one. 


Figures A3 and A4 are recordings that show the ability of the 
thermocouple system to produce the identical results from the same inputs 


on different runs at different times. 
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APPENDIX B: STANFORD UNIVERSITY RADIANT HEATING STRUCTURAL TEST FACILITY 


I. INTRODUCTION 


The need for aerodynamic heat simulation in structural testing of 
supersonic aircraft is well recognized. Thus, the authors undertook 
the design and construction of a four channel radiant heating test facil- 
ity at Stanford University during the summer of 1962. The system was 
designed to meet the following requirements: a high intensity, rapid 
response heat source; an adequate reflector system to direct the heat; 
an absorptive finish on the test specimen to permit maximum utilization 
of available heating power; a practical technique for surface temperature 
determination; a system of radiation measurement; and a computer control 


system to correctly vary the heat input to the test specimen. 


II. SYSTEM LOGIC 


Figure B-l is a block diagram of the control svouen! for the radiant 
heating test facility. Voltages analogous to h and ae are generated 
and then utilized by the computer along with a voltage analogous to ve 
for solution of the Newtonian convective heat transfer equation, 

q = h(t - t.) . Another voltage (a) representing radiation fron 
external sources (sun, atomic blast, etc.) is generated and then added 
to q ina summing circuit to give the desired heat flow signal, qa: 
The comparator output, c, is the difference between qa and ae 

the measured heat flow. The control signal causes the power regulating 
device to vary the power supplied to the heaters. The heat flow measure- 
ment system closes the control loop, causing the correct heat input to 
the test structure to be maintained. The measurement of heat input to 
the test specimen is made by either of two systems. The operational 7 
details of these devices are explained briefly in Section III and in 
detail in References [1] and [2]. 


III. SYSTEM DETAILS 


4 


Certain requirements are essential in a practical feedback control 
system for the simulation of aerodynamic heating. The requirements set 
for this facility are to: 
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(1) accurately compute desired heat flows for simulation of speeds 


up to Mach 6. 
(2) be capable of measuring the heat inputs to the test specimen. 
(3) be dependable. 
(4) be safe to operate. 
(5) be inexpensive as nosei neem 


(6) be simple to operate and program. 


The operations of the function generator and computer shown in 
Fig. B-l are combined in the digital to analogue computer control systen. 
This unit consists of a punched tape reader, an operational amplifier 
and resistance bank, two peak followers and associated capacitors, a 
multi-bank stepping switch, a summer, and a chopper stabilized amplifier. 
(Fig. B-2) The details of the individual components are given in Ref- 


erence [1] and only the modus operandi is discussed here. 


Values of h and ce are obtained from the Stefan-Boltzmann law. 
A binary coded program of these time dependent functions is made on a 
punched tape. This tape is placed on a reader and used in the time 
Shared system as follows: an analogue voltage proportional to q.. is 
generated from a resistance network and a reference voltage in conjunction 
with an operational amplifier. This signal is temporarily stored. The 
next command is Lc which is treated in the same way. The output of 
the temperature transducers, te » is amplified by a chopper stabilized 
amplifier and then subtracted from the stored eer analogue in a summer. 
The voltage difference ae zs t.) is used as the reference voltage for 
the next stage of the operation. This consists of giving the command, 
read h , and by virtue of the fact that the reference voltage is 
analogous to (t - t_) the output of the system is analogous to the 
product n(t_ - t.) - Simultaneously, q and q_ are added in a 
Summer and this voltage, analogous to qq » is stored in a capacitor. 
These steps are repeated for channels two, three, and four. The digital 
to analog computer evaluates a new value of qa for each channel every 
1.2 seconds. Thus, by time sharing, by digital to analog conversion, 
and by unique. computer pec iiones Wc @ simple and accurate means for 


obtaining a qa history is realized. 
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There are two systems available for closing the control loop. Both 
systems meaSure the rate of heat flow to the test specimen and therefore 
may be used interchangeably. One system, the classic system, consists of 
special transducers. These essentially consist of thermocouples in con- 
junction with substantial heat sinks. Ina radiant field the thermo- 
couples reach an equilibrium temperature proportional to the radiant 
intensity. Each transducer is two thermocouples connected so that their 
signals oppose. One is subjected to the radiation from the heaters and 
the other to the reflected radiation. The output is therefore proportional 
to the difference between the incident and reflected radiation. This is 
qn exclusive of convective losses. The heat sinks for the transducers 
are cooled metallic blocks. The instruments are attached directly to the 


Sides of the heater arrays. They are calibrated in a calorimetric manner. 


The other system, one especially devised for this project, consists 
of eighteen thermocouples arranged to obtain temperature profiles in the 
test specimen. A high speed stepping Switch sequentially connects the 
thermocouple signals to an analogue computer which solves the appropriate 
heat flow equation. The computer is time shared by the 4 stations. The 
latter system is more versatile than the former since it measures net 
heat flow irrespective of the nature of the source. Thus convective losses 
in simulation are accounted for. An additional refinement to the computer 
permits compensation for radiant losses which the structure would normally 
experience in flight. Figure B-3 is a photograph of the feedback con- 


puter system. 


Voltages in the control system are scaled such that at any time the 
prescribed de level is twice the anticipated qn level. <A summing 
circuit subtracts ote from qa and the output (c) is used to control 
the power regulating devices. Due to the two-to-one scaling the control 
Signal consists of a set-point voltage plus an error voltage. Thus, the 
heaters are controlled in such a manner that the deSired and actual heat 
inputs are equalized. This occurs when the error voltage level converges 


to zero and the control signal converges to the Set-point voltage. 


Main power is supplied through two 300 KVA, 480 volt, 3 phase, 
General Electric transformers. The power regulating devices consist of 


75 KVA, 480 volt, 3 phase General Electric saturable core reactors and 


= 3l i as 


| 


— > ae _— iin 








their associated control units. Transient operation of the power equip- 
ment on an overload duty cycle makes possible peak outputs of two to 
three times rated value. Heating capacity of more than one megawatt is 
therefore feasible. Figure B-4 is a photo of the saturable core reactor 


Wnts. 


Tne basic heating elements are General Electric T.5 infra-red lamps. 
(See Fig. B-5 ) They consist of a tungsten filament supported by tanta- 
lum disc spacers in a sealed quartz tube 3/8 of an inch in diameter and 
12 inches in length. Each unit is rated at 1000 watts at 220 volts. 

They are capable of withstanding the 480 volt peak output of the power 


supply and under these conditions develop approximately 3 kilowatts. 


The heaters are arranged in groups of six. Figure B-6 is a photo- 
graph of one such unit. The heater arrays are connected in delta to the 
three phase power from the saturable core reactors. Each unit covers 
about one third of a square foot. Thus, the maximum radiant intensity 
possible with the present arrangement is about 60 kilowatts per square 
foot. Two thousand watt heaters of the same external size are avail- 
able when higher intensities are required. Peak output of these is 6000 
watts per heater. The heater arrays are quite adaptable and may be 
bolted to any suitable frame to surround a test structure. Special 
arrangements may easily be devised for special purposes, such as the 


heating of curved surfaces. 


Figure B-/ shows an overall view of the test facility. The main 
items of equipment and the channels of information flow are outlined on 
the photograph. All control, computing and recording equipment is located 
in the control room. This room is on the second floor of the laboratory 
Overlooking the test site. Large windows provide for visual observation 
of the structure under test. The remote location of the control room 
was dictated by safety considerations and provides a convenient, compact 


center for operation of the entire facility. 


IV. CONCLUDING REMARKS ‘ 


The design, construction, and operation of this test facility has 
proved that the concept of component time-sharing is practical in certain 


types of control systems. It is concluded that such time-sharing is an 








economical proposition in which savings in equipment costs are balanced 
against slightly reduced system performance. The facility is such that 
it can be readily installed in any small research laboratory where cost 


is of prime importance. 
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Saturable Core Reactors 
Fig. B-4 


Heater Array 
Fig. B-6 
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Infra-Red Heater 
Fig. 
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Overall View of Test Facility 
Fig. B-7 
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Control Room 
Fig. B-8 
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